Suppression of Hyperfme Dephasing by Spatial Exchange of Double Quantum Dots 



(N 

o 

< 

ov 

(N 



c3 



T3 

O 
o 



> 

in 
o 

od 
O 
(N 



X 




FIG. 1. Suggested electrode geometry for a rotating double- 
quantum-dot qubit with top and bottom gates in different 
shades of gray. Exchange gates via real space rotation, as 
opposed to tunneling, are expected to strongly reduce the 
qubit sensitivity to charge noise. 
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We examine the logical qubit system of a pair of electron spins in double quantum dots. Each 
electron experiences a different hyperfme interaction with the local nuclei of the lattice, leading 
to a relative phase difference, and thus decoherence. Methods such as nuclei polarization, state 
narrowing, and spin-echo pulses have been proposed to delay decoherence. Instead we propose to 
suppress hyperfme dephasing by adiabatic rotation of the dots in real space, leading to the same 
average hyperfme interaction. We show that the additional effects due to the motion in the presence 
of spin-orbit coupling are still smaller than the hyperfme interaction, and result in an infidelity below 
10~ 4 after ten decoupling cycles. We discuss a possible experimental setup and physical constraints 
for this proposal. 



I. INTRODUCTION 

In recent years there has been great interest in the 
prospect of using scalable solid state devices to imple- 
ment quantum two-level systems (qubits) for potential ap- 
plications such as quantum computation. One promising 
candidate for a qubit is a pair of electron spins in quan- 
tum dots, which forms a fault-tolerant subspace that is 
immune to collective decoherence .|j However, each elec- 
tron is still subject to the local hyperfme interaction from 
the nuclear spins of the lattice, which leads to dephasing 
of the individual electron spinsB There have been sev- 
eral proposals to suppress this dephasing such as nuclear 
polarization. jd state narrowing,El and spin-echo pulse 
correction.Qi3 While improved coherence has been exper- 
imentally demonstrated using these techniques, the co- 
herence times desired for applications have proven very 
difficult to achieve. For example, pumping methods have 
been used to partially polarize nuclei, but the. nearly full 
polarization needed has yet to be achieved. Bulij 

Promising results have been shown using spin-echo se- 
quences pfchrough the exchange interaction between two 
spinsEMLj The exchange is controlled by lowering the 
tunneling barrier between the two quantum dots using 
quickly-controlled electric gates. This leads to Rabi os- 
cillations; a single 7r-pulse corresponds to exchanging the 
two spins. Sequences of such pulses effectively couple 
both spins to the same average hyperfme interaction re- 
sulting in improved coherence times. While several echo 
sequences can be performed using this exchange, it is 
likely that the electrical gate noise and spatial variations 
in the Ovexhauser field remain the dominant sources of 
dephasing.Ej 

Rather than relying on interdot tunneling, we propose 
using a spatial exchange of the two quantum dots, allow- 
ing the two electrons to traverse the same path, spend- 
ing the same time coupled to the local nuclei, as shown 
in Figure ^. Compared to exchange via tunneling, ide- 
ally, this should eliminate the effect of the electrical gate 
noise. On the other hand, in the presence of spin-orbit 
coupling, such a motion of electrons may introduce ad- 
ditional errors. In this manuscript we analyze how these 
errors depend on the parameters of the motion and dis- 



4k 

_ 371 

a 

a 

% 27C 

<2 



Single Direction/ 



Alternating \ 



0.5 1 1.5 2 

t/T 

FIG. 2. Time dependent adiabatic trajectories used in 
the simulations. Plotted is the position of the first dot 
parametrized in terms of the angle 9 as a function of time 
t. The angle-dependent positions were defined as a sum of 
properly scaled and shifted hyperbolic tangents. Single direc- 
tion rotations suppress the effect of a static Overhauser field 
but not of a time- varying one. Longer rotation sequences like 
alternating forward-and-back suppress the effect of a linear in 
time Overhauser field. 



cuss the constraints and possible parameters for potential 
implementation of this proposal. 

We note that a similar setup with spatial ex- 
change of electrostatically-defined quantum dots has 
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been discussed .in relation to holonomic quantum 
computation. oEj However, the corresponding coherence 
estimates have been done in the absence of a magnetic 
field. Here we focus on the parameter range characteris- 
tic of double-quantum-dot qubits and take into account 
typical magnetic fields of order > 0.1 T. 

The outline of the paper is as follows. We define the 
Hamiltonian of the double-quantum-dot qubit with spa- 
tial exchange in Sec. |n|, derive the effective spin-only 
Hamiltonain for a single spin in a moving quantum dot in 
Sec. HI, and calculate the single-qubit fidelity associated 
with a sequence of double-dot rotations in Sec. [IVL We 
discuss simulations of a possible protocol in Sec. [V| and 
the constraints and corresponding characteristic time and 
distance scales in Sec. VI, and conclude with Sec. VII. 



II. SETUP 

We envision a qubit formed by a pair of quantum dots 
electrostatically defined in a III-V semiconductor (e.g., 
GaAs/AlGaAs) heterostructure using a system of top 
and bottom gates similar to that illustrated in Fig. [|. 
We take the parameters of the dots to be similar to the 
experimental ones. Ejllj Specifically, each dot contains a 
single electron, with a typical dot-size quantization en- 
ergy fkud ~ lmeV. The qubit is defined as the subspace 
of the singlet and m s = triplet states of the two elec- 
trons. The triplet degeneracy is removed by a uniform, 
constant magnetic field Bq of at least ~ 0.1T applied per- 
pendicular to the sample, which creates a Zeeman gap of 
A ~ 2.5 ueV. The electrons interact with N ~ 10 6 spins 
of the lattice nuclei, leading to a different local hyper- 
fine interaction for each electron. We will approximate 
this as a Zeeman interaction with a random, fluctuating, 
non-uniform magnetic "Overhauser" field ~ lmT. 

To prevent dephasing, both quantum dots will be 
moved along the same trajectory in a time much shorter 
than the relaxation time of the nuclear spins, t nuc ~ 
100 us, so we can assume that the Overhauser field is 
quasi-static. In order to reduce the sensitivity to charge 
noise the distance between the dots must be significantly 
greater than the size of each quantum dot, a ~ lOOnm. 
Due to this spatial separation between the dots, we can 
treat the Hamiltonian of each electron separately, 

H = H d (r Q ) + H Z + V z (r,t) + H so + V(r,t). (1) 

Here the dot Hamiltonian is given by 



#d(ro) = |-+t/(r-ro(i)), 
2m 



(2) 



with the canonical momentum p = P + eA/c and the 
confining potential U centered at r = r (t); the Zeeman 
Hamiltonians for the externally-applied and Overhauser 



fields are respectively 

nj, — — — rs • rr , 
V z (r,t) = ^-B N (r,t)-<r, 
and the spin-orbit Hamiltonian is given by 



H. 



so 



(3) 
(4) 

(5) 



The last term, V(r,t), accounts for additional effects 
originating from disorder, variation of the dot poten- 
tial as it moves due to imperfections of the confining 
potential, as well as phonons. The spin orbit coupling 
coefficients C,-j in Eq. (||) incorporate both Dresselhaus 
(originating from the lack of the inversion symmetry of 
the lattice) with C vv = —C xx = /3, and Rashba terms 
(structural inversion asymmetry due to the quantum 
well) with C xy = —C yx = a. This specific form as- 
sumes that the growth of the semiconductor heterostruc- 
ture and t h p Q ua ntum well asymmetry are in the positive 
z direction,E30 so all the matrix elements that involve z 
are zero. 

For numerical estimates we will use the effective elec- 
tron mass to ~ 0.067m c , and assume a ~ /3 with values 
ranging from 10 3 to 10 4 ni/s, as appropriate for typi- 
cal GaAs heterostructuresJiS Based on the above val- 
ues, we have mf3 2 <C Q/J-bB <C hu>d, so we will ignore 
terms quadratic in the spin-orbit coupling for the follow- 
ing analysis. 



III. EFFECTIVE SINGLE-DOT HAMILTONIAN 

We now go into the moving reference frame of the dot 
using the translation operator 

*(<) -> e-* p ' r °W*(t), (6) 
H -> e i p - r °( t ) J ffe^ p - r ° (t) - v (t) • P, (7) 

which introduces the additional term proportional to the 
dot's velocity, vq = i"o, and removes the i"o from the con- 
fining potential, f/(r — ro) — > U(r). This also affects the 
vector potential, A(r) — > A(r -I- r ), in the dot Hamilto- 
nian and spin-orbit term, which can be reversed with an 
appropriate gauge transformation. We select the sym- 
metric gauge, A(r + ro) = |B x (r + ro), and transform 
A — >• A + V/ with / = — ir • (B x r ) which results in 



* ->• * cxp(- — /), 
nc 

P -> P -V/. 

c 



(8) 
(9) 



These transformations introduce two additional terms 
in the time-dependent Schrodinger equation ihd^S/dt = 
i/vf . The first term, — (e/2c)vo • (B x ro), arises when 
Eq. (||) is substituted into the — v (t) • P term in Eq. (Q). 
The second term, — (e/2c)r ■ (B x vq), appears on the left 
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hand side as a result of taking the time derivative of the 
exponent in Eq. (|J). These two terms can be moved to 
one side of the equation and combined using the cyclic 
property of the mixed product: 



- — v • (B x r ) 

2c 



6 6 

— r • (B x v ) = -tt v o • B x (r + r ) 

AC AC 

(10) 

which is precisely the vector potential term in — v - p. 
This leads to the moving-frame Hamiltonian 

H = H d + H z + H so + V(t + r (f), t) 

+ Vz(r + r (t),t) -v -p. (11) 

Following Loss et alEl, we now perform a canonical trans- 
formation H -> e s He~ s ~ (1 + S)H(1-S) = H+[S,H], 
where S is anti-Hermitian and chosen to eliminate the 
original spin-orbit term. We split S = So + S% such that 



[So, Hd] + Hso — 0, 
[Si, + [So, Hz] = 0, 



and choose 



So 



(12) 
(13) 



(14) 



which satisfies Eq. (12). This can be verified, noting that 



Sq has no momentum dependence so it clearly commutes 
with the confining potential, and 



%m isyii 4 V 

—a % C n r 3 , — 
n 2m 



-a % C iJ [r 3 ,p k ]p k 



= -o l C lJ p>. 
With So known, we use Eq. (O) to define Si, 



(15) 



[Ha, Si] — [S , H z ] 



2h 



-C J r J [a , a \B 



^(B„xQ).,, 



(16) 



where Q % = {2m/K)C il rK This equation can be written 
in terms of the electron's orbital states \n) in the dot 
potential, 

(n\[H d ,Si]\m) = (Si) nm (E n - E m ) 

= 9 -^-a ■ (B x (Q) nm ). (17) 

As long as the relevant dot quantization energies are non- 
degenerate, E n ^ E m , we have 



to \ _ 9VB (B x (Q) nm ) 



w 

vv ri 



(18) 



where we defined 



Wnr. 



Bp x (Q) r 

E n — E n 



(19) 



Expanding the canonical transformation to first order in 
the spin-orbit parameter contained in Sp and Si , we ob- 
tain the transformed Hamiltonian 

H~H+ [So, V z ] - [So, v • p] + [S U V] 

+ [Si,H z ]- [Si , v ■ p] + [Si , V Z ] (20) 

Using the definition of So, the first two commutators are 



img/iB 
2h 



C ij r i [a i ,a k ]B^ = - 



mgnB 



C t3 r 3 (e M B&o 1 ) 



5MB 



(B N x Q) ■ er 



(21) 



and 



-^C ij [r j ,p k H = ma^vi = ±Q V ■ a, (22) 

respectively, with Q\ = 2mC ZJ v 3 Q defined analogously to 
Q l . We now define the effective spin Hamiltonian by pro- 
jecting onto the orbital ground state, Hs = (0|-ff|0). This 
allows us to express the remaining commutators in the 
transformed Hamiltonian using the definition of Si . The 
first commutator involving Si in Eq. ( p0| ) is simplified 
by explicitly writing out the commutator and inserting a 
complete set of states, 



(0|SiF|0)-(0|FSi|0) 

= J2(0\Si\n)(n\V\0) 



n>0 



)|y|n)(n|Si|0) 

= /J(Sl)onKlO ~ Vo,i(Si)„o- (23) 
n>0 

We assume V to have no momentum dependence, so it 
commutes with Si and we can reverse the order of the 
second term above. Since Si is anti-Hermitian, while V 
is Hermitian, 

<0|SiV|0)-(0|VSi|0> 

—2 y^(Si)o TO V^,o 

n>0 



(24) 



n>0 



This technique can be used on the remaining terms in 
Eq. (p0|). The second term involving Si in Eq. (20) van- 
ishes because 



(H Z )n0 



tr- (n|B o |0> = 0. 



(25) 



However, the final commutator in Eq. (|2fj) contains the 
terms 



(Sl)p„(Vz)„p — (Vz)p n (Sl)„p, 



(26) 
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which do not commute because they contain two spin 
terms, but can be treated using the spin identity 

O • a)(er • b) = a- b + ier ■ (a x b). (27) 

The first term in Eq. (|2^) becomes 

(^p) 2 {Won • (B N ) n0 + i<r ■ (Won x (B N )„ )}, (28) 

and the second term looks quite similar, except the anti- 
commutator of the cross product causes the spin depen- 
dent term to cancel with the one above, while the spin- 
indepedent terms is doubled. The Hamiltonian contains 
several of these spin-independent terms that can be taken 
as constants. The effective spin Hamiltonian, up to a 
constant, can now be written simply as 



where 



, , - 9flB "R 



(29) 



(30) 



is the Larmor frequency, and the time-dependent term 
u!i = u!i(t) is 



h 

1 

2fi 



Bn + (Bn)ok x (Q)„o] 



4.9Mb 



WonKo 



Tl>0 



1 Q v - ^2 W 0n( v • Pno)- 

n>0 



(31) 



If we ignore the phonons and approximate the Over- 
hauser fields as static, the time dependence of the terms 
in the first line of Eq. ( |3~i| ) comes only from the posi- 
tion, parameterized by the known trajectory of the dot, 
r (t) . Similarly, the time-dependence of the terms in the 
second line comes from both the position ro(t) and the 
dot velocity v (t) [in fact, these terms are all linear in 
components of Vo(i)]. The simple spin Hamiltonian ( p9| ) 
is the key result of this derivation; it is correct to linear 
order in the spin-orbit couplings. 



IV. AVERAGE FIDELITY 

A. General expression 

In order to analyze the implications of the additional 
terms in the effective spin Hamiltonian (p9|), we need to 
take into account that the qubit is actually formed by two 
electron spins. It will be convenient to assume that the 
dots' velocities are small compared to the speed of sound 
»o C s ~ 5 X 10 3 m/s. Then the phonon effects should 
decouple from the dots' motion and can be approximated 
as contributing to the same "intrinsic" decoherence times 
T\ , T2 as one would have without the motion. The effect 
of such decoherence terms on dynamical decoupling has 



been considered in detail in Ref. |l9| in the following we 
assume that these decoherence times are large compared 
to the characteristic period T of the dots' motion and 
therefore can be ignored. 

In the absence of phonons, and approximating the 
Overhauser field as classical, the time evolution of the 
two-spin wavefunction with N — 4 components can be 
characterized by the unitary matrix U (t) . The qubit 
subspace Q is formed by the m z = component of the 
two-spin wavefunction; it has M — 2 dimensions. The 
standard assumption is that, at the beginning of the ex- 
periment, the two spins are initiated in a pure state 
which belongs to the qubit subspace, \if}) € Q. Therefore, 
when computing the average fidelity, we need to average 
only over the original wavef unctions in Q. 

More generally, consider an M-dimensional subspace Q 
of an iV-dimensional Hilbert space H. Let us introduce 
an N x M matrix T whose columns are formed by the 
components of orthonormal vectors forming a basis of 
Q. Then, the components of an arbitrary wavefunction 
\ip) E Q are a linear combination of the columns of T; 
namely ip — Tip, where tp is an M-dimensional column 
vector, \\ip\\ = 1. The corresponding density matrix can 
be written in this basis as po = Tipip^Ti. The fidelity 
corresponding to the evolution matrix U is 

F = Tr{ Po Up tf) = (<piW<p)(<pW<p), (32) 

where W — T^UT can be thought of as the projection of 
U onto the subspace Q. The average fidelity in the sub- 
space can now be calculated using the averaging identities 
for components (pi of the normalized wavefunction \p) 



{<Pi<p*} = Sij/M, 



(iPiP*ip k ip^) = 



SijSki + SuSjk 
M 2 + M 



This leads to the average fidelity 

\TtW\ 2 + Ti(WW^) 



(F} = 



M 2 + M 



(33) 
(34) 

(35) 



For the special case of the qubit formed by the singlet 
and m = triplet states of two spins, assuming no inter- 
dot tunneling, we can take the net evolution matrix as 
the Kronecker product of evolution matrices correspond- 
ing to the two qubits, U = U\ ® Ui- Further, it will be 
convenient to decompose each single-spin matrix in the 
interaction representation with respect to the precession 
in the net effective magnetic field along the z-axis, 



Ui — UoiSi, Uoi 



where 



<Pi(t)=w t+ / dt'uj z u (t'), 



(36) 



(37) 



ujq is the Larmor frequency, wfj(t), i= 1,2 [cf. Eq. (pfi|)] 
are the two dot's effective perturbing fields in the z- 
direction, and the matrices 

Si = e-^i-itfv-/^ 1 = 1,2, (38) 
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are parametrized as rotations by angle (pi = \<pi\ around 
the unit vectors <pi, with extra phases ji. These rota- 
tions come entirely from transverse, fi — x,y, compo- 
nents of u)^ in the rotating frame, largely due to the 
Larmor frequency. Since the Larmor frequency is large, 
the additional rotation angles are expected to be small; 
we expand the average fidelity ( ^5|) to quadratic order in 
components of cp i , 



(F) = l- f - h - f 2 - / 2 J 



with the infidelity terms 
2 

3 
1 

3 

ft= 2 - 



fo 
h 



sin 2 (A^/2), 
(<f> z 2 - 4>\) sin(A^), 

0l) 2 COs(A( y 5), 



cos(A<,5) 
~T2 



±\2i 



utr + (#) 



(39) 

(40) 
(41) 
(42) 
(43) 



Note that, as expected, the fidelity only depends on the 
differences Aip = tpi—tp\ and <p 2 — 4>f of the two precession 
angles around the z-axis. To the same quadratic accuracy 
in the small angles, we can also write 

fo + fi + /I = \ [1 - cos(A^ + <&-<&)l (44) 

which only depends on the total phase difference, and is 
exact in the limit <&r =0. 



B. Rotating- frame approximation 

We now return to the analysis of the single-spin Hamil- 
tonian (p^). The Larmor frequency, u> > 4 x 10 9 rad/s, 
is the dominant term, that is, wi <C wo- Given that we 
have excluded the phonons, we can also assume that the 
dot trajectory is such that the time dependence in oj\(t) 
is slow on the scale of cjo- This implies that we can use 
the interaction representation Eq. (p6[), where S(t) is the 
slow part of the unitary evolution operator; it obeys the 
equation 

ihS = H int (t)S, 5(0) = 1, (45) 
where we temporarily omit the dot index, and 

-H^t) = |«i(t) • U$(t)trU (t) 

= ^{a x [u*{t) cos (p(t) +u)\(t) sin ip(t)\ 

+a y [u\ (t) cos (p(t) - uj* (t) sin ip(t)} } (46) 

is the perturbing Hamiltonian in the interaction repre- 
sentation. 

The remaining analysis is performed perturbatively 
with the Magnus (cumulant) expansion, 



where the first two cumulants are 
C (1) (i) =-\ I «ftifitat(ti), 



6 2 \t) 



1 

2fr? 







dti dt 2 [H int (t 2 ] 

0<t!<t 2 <t 



(48) 

,^int(ii)]. (49) 



We perform the integration explicitly to leading order in 
oji/ljo, also assuming that the time-dependence of lo\ is 
slow on the scale of wq: 



{*) = —{(J x {a sin <p(t) + b Q - b cos ip(t)} 
- a + acosip(t)]X, 



2w 
-a y [b sin ip(t) 



1 /"* 

iCW(t) = a z / dh [a< 

4w Jo 



(t^+b 2 ^)). 



(50) 
(51) 



Here we denoted a = cuf (t), b = u)\(£), and oq, the 
corresponding values at t = [</?(0) = by definition]. 
Eq. (p9|) gives the expression for the qubit infidelity 



1 - (F) = f z + ft + f 2 



r 

n 



= - [1 - cos (A^ + 8 2 -5 x ) 



(52) 
(53) 



■ COS(Aly9) 

12cj 2 



[{Bi - b 0i f + (Ai - a 0l ) 2 ] , (54) 



where the index i = 1, 2 refers to the two spins, Ai = 
ai cos tp(t) + bi sin ip(t) and Bi = bi cos ip(t) — a,i sin ip(t) 
are the rotated components of the transverse angular- 
velocity vectors (u)f,uif) for the corresponding spins, and 
the additional phases Si are given by the integrals. 



5i = 



dt' 



[<4(f')l 



(55) 



One immediately recognizes the additional phases in 
Eqs. ( |53"| ) and (55) as the effect of level repulsion, or 
equivalent ly as the gap for the spins, driven adiabatically 
by the total instantaneous magnetic field oc [(wo+w^) 2 + 
(u;^) 2 ] 1 / 2 . Then, Eq. (|||) can be interpreted as the ef- 
fect of the basis change between the original quantization 
z-axis and the direction of the instantaneous magnetic 
field; it is similar in nature to the initial decoherence as- 
sociated with any periodic decoupling sequence, see, e.g. 
Ref. |(| 

The terms of the next order in 1/ujq expansion, omit- 
ted in Eqs. (|50| ) and (pT|) , include the trivial correc- 
tion to Eq. ( |51| ) oc J dtu>l[uji] 2 /coq, as well as the ge- 
ometrical phase oc J dt' W[u)f , oj v ], where W[x, y] = 
W[x{t), y{t)} = x{t)y'{t) - x'\t)y(t) is the Wronskian. 

Since Aip = f Q dt' [u>* 2 (t') ~ ^u(t')], the term f z does 
not contain any rapid oscillations at the Larmor preces- 
sion frequency ujq, while the terms can be averaged 
over the period of Larmor precession by replacing A 2 +Bf 
with {oJii) 2 and dropping all of the terms linear in Ai, Bi, 



(47) 



7t = 



2 + cos(0| - <pl) 



3w 2 



{UM 2 + UM} 2 }- (56) 
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C. Sequences 

Overall, the dynamical decoupling should be designed 
to null the difference between the accumulated phases 
<Pi + Si of the two spins which suppresses the main contri- 
bution to infidelity, see Eq. ( f53|) . For a static Overhauser 
field, this can be achieved just by ensuring that each spin 
spends the same amount of time at each position, e.g., 
via the solid adiabatic trajectory in Fig. M. This is also 
sufficient to suppress the effect of the velocity-dependent 
terms in the second line of Eq. (Q) . A more complicated 
set of dot rotation involving motion in both directions, 
e.g., see the dashed trajectory of Fig. are required to 
suppress a time-dependent Overhauser field. 

To analyze the effect of a sequence of 7r-rotations in the 
presence of a time and position-dependent Overhauser 
field, consider its Fourier expansion at a position on the 
trajectory parametrized by the rotation angle 9, 

B z (9,t) =A (t)+y'A m (t)cosm8+B m (t)smm6. (57) 



Only the difference between the fields corresponding to 
the two dots (located at 9 and 8 + n) is relevant for the 
infidelity Eq. ([33]). This leaves only the terms with m 
odd in the Fourier expansion (|57|). 

For a term with cos m9 (an even function of 9) a se- 
quence of 7r rotations acts the same way, independent of 
the direction. It is easy to check that with an equidis- 
tant sequence of rotations centered at Tq/2, 3Tq/2, . . . , 
(2s + l)Xb/2, where the number of rotations s is even; 
any time-independent and linear in t contributions to 
A m (t) are suppressed, but a quadratic term would gen- 
erally remain. Unlike the usual dynamical decoupling 
problemQ'Eil, it is not generally possible to suppress the 
quadratic term of A m (t). 

The rotation direction starts to matter for a term with 
s'mm9 which is an odd function of 9. Here a sequence of 
7r rotations in the same direction picks up a sum of contri- 
butions from consecutive time intervals with alternating 
signs, suppressing the time-independent contribution to 
B m {t) but not the linear contribution. As an alternative 
prescription, a symmetrized sequence of two forward ro- 
tations by angle 7r, followed by two rotations in the oppo- 
site direction can be used to suppress the effect of the lin- 
ear in t term in B m (t). Generally, it is possible to design, 
more complicated sequences analogous to concatenated!!!! 
or Uhrig'su sequences to suppress the effect of any fixed- 
degree polynomial in time B m (t). However, we do not 
expect this to be useful since the quadratic time contri- 
bution of A m (t) would still remain. 



V. SIMULATIONS 

We corroborate our analytical results by simulating the 
two-spin unitary evolution with the effective Hamilto- 



by the rotation angle 9 = 9i(i)\ the other dot is assumed 
to have the symmetric position, 9 2 {t) = 9\{t) + w [see 
Fig. U for samples of actual trajectories.] The position- 
dependent terms in the first line of Eq. (|3l| ) are simulated 
in terms of a three-component correlated magnetic field 
B(#) drawn from the Gaussian distribution with zero 
average and the correlation function (B li (9)B lJ {9')) = 
o'f l 5 lll/ '&(9 — 9') (no implicit summation in /i, v = x, y, z), 
where 



-{9-2-nm) 2 /I 2 



nian (29). Specifically, we parametrize the dot trajectory 



is an infinite sum of Gaussian functions (which can also 
be represented in terms of the Jacobi theta function). 
These are obtained by applying a Gaussian filter to a 
discrete set of uncorrelated random numbers drawn from 
the Gaussian distribution, and using the standard cu- 
bic spline interpolation with the result. To simulate the 
components of time-dependent magnetic field B(r,t), we 
used explicit order-r spline interpolation between several 
such angle-dependent functions, where r = 1,2. 

For all simulations, we chose the time units correspond- 
ing to the Larmor precession period To = 2tt/ujq and the 
correlation time of the Overhauser field 4- 10 4 to with each 
component of its r.m.s. value corresponding to rotation 
frequency | 2 ) = 0.025/tq. The adiabatic trajectories 
of the dots were simulated using a sum of appropriately 
shifted hyperbolic tangents, scaled so that the dot is in 
motion during approximately half of the protocol. The 
leading velocity-dependent term in Eq. ( |3l| ) was simu- 
lated using the corresponding derivatives and the param- 
eter Q v /h — 0.075/to, assuming equal contributions from 
the Rashba and Dresselhaus parameters. 

For the case of the static Overhauser field, see Fig. 0, 
we see that the average infidelity is dominated by the 
contribution from the z-component of the field, and that 
the infidelity nearly vanishes at the end of the spatial ex- 
change protocol. For a linearly time-interpolated Over- 
hauser field, the infidelity increases over several cycles 
of the single-direction tt pulses, Fig. [|a, but maintains a 
low value ~ 1CU 5 after alternating between a sequence of 
two forward rotations, followed by two rotations in the 
opposite direction, Fig. J^Jd. However, for the quadratic 
interpolation, Fig. |p, the infidelity gradually increases 
even for the alternating protocol, though it stays below 
10~ 4 for several cycles. 



VI. POSSIBLE EXPERIMENTAL SETUP 

Aspects of our proposal, such as the precise construc- 
tion of few-electron quantum dots in III-V semiconduc- 
tor heterostm<j^ps_have already been demonstrated in 
experiments. liSlO't^lLj However, the precise adiabatic ro- 
tation required in our proposal may be quite difficult to 
accomplish experimentally. We now discuss possible de- 
sign implementations for an experimental realization, as 
well as physical constraints. 



7 




399.5 400 400.5 

t/T 

FIG. 3. Simulated qubit infidelity 1 - (F) [Eq. @] in the 
vicinity of the first full rotation period of the double-dot qubit 
at t = T. Position-dependent magnetic field B^iO) is as- 
sumed static, and the rotation period T is chosen commensu- 
rate with the Larmor frequency, uioT/2-n — 400. Dashed line: 
only B z {9) is included; the infidelity minimum is exactly at 
t = T, in agreement with Eq. which is exact in this sit- 
uation. Dotted line: only the transverse components 5^(6), 
fi = x,y are included. The infidelity minimum is slightly off 
the commensurate time t = T due to the terms not included 
in Eq. (M). All three components of the field B^(6) are in- 
cluded for the red solid line. 



gular direction could be accomplished by placing appro- 
priately chosen time-dependent voltages on the "wedge- 
gates" on both sides of the electron (Fig. [l]). Several of 
these wedges will be needed to accomplish the smooth 
and adiabatic trajectory needed. The combined use of 
wedge and circular gates may require gates on both the 
top and bottom of the sample. The typical confining 
potential is approximated by U ~ muj^a 2 , which only re- 
quires reasonable gate voltages on the order of 100 mV. 

The averaging of the hyperfine interaction is only valid 
in the quasi-static approximation of the Overhauser field, 
which requires that T -C t nuc ~ 10~ 4 s. This places a 
lower bound on the velocity, while there is also an upper 
bound necessary to ensure that the dot rotation does not 
produce a Lorentz force that causes the electron trajec- 
tory to deviate by more than a correlation length. This 
results in the restriction, 10 _1 <C v <C 10 5 in m/s. In 
our estimates we also assumed v to be small compared 
to the speed of sound, s~5x 10 3 m/s. With a/ro ~ 15 
and T ~ 1 /is, v ~ 10 m/s easily satisfies this condi- 
tion and, according to our simulations, should result in 
the infidelities lower than 10 -4 for significant time-scales. 
The velocity could potentially be increased to allow more 
operations to be performed before the states decohere. 
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10 15 20 25 30 35 
t/T 

FIG. 4. Simulation results for fidelity measured at the end 
of each cycle, for a linear time-dependence of B^ with (a) 
forward, (b) forward-back dot rotations [cf. Fig. B, as well as 
(c) forward-back for quadratic time-dependence. 



As discussed above, our proposal is suitable in ma- 
terials with relatively weak spin-orbit coupling such as 
GaAs/AlGaAs heterostructures. We believe the elec- 
trons in a 2DEG could be confined in the radial direction 
by creating a circular depletion layer by placing electro- 
static gates with a constant voltage in the center and 
outer edge of the circle as sketched in Figure [|. Since 
we wish to suppress the tunneling, the normal inter- 
dot spacing should be much larger than the dot size, a, 
so we assume that the circular trajectory has a radius 
ro ~ 15a ~ 1 um. Confinement and rotation in the an- 



VII. CONCLUSION 

We analyzed the real-space exchange of quantum dots 
as a possible substitute for the tunneling exchange. Ide- 
ally, exchange eliminates the hyperfine dephasing from 
the Overhauser field parallel to the applied field, leav- 
ing only the smaller effects from the in-plane field. The 
real-space exchange accomplishes the same suppression 
of the hyperfine interaction, but avoids the problematic 
sensitivity to charge noise present in exchange via tun- 
neling. While spatial exchange does introduce additional 
effects such as spin-orbit coupling, simple tricks like using 
pairs of 7r rotations in alternating directions can be used 
to suppress these so that the decoherence is still dom- 
inated by the hyperfine interaction. In particular, this 
field only enters as the small ratio of the average in-plane 
Overhauser field to the externally applied field. Perhaps 
the simplest way to suppress the hyperfine interaction in 
this approach is to reduce this ratio by increasing the 
externally applied field. 

In addition, this spatial exchange is compatible with 
some of the methods already being attempted such as 
nuclear polarization via pumping. If the hyperfine inter- 
action can be further suppressed, the next largest contri- 
bution from our spatial exchange approach would come 
from the disorder of the sample or the electron-phonon 
coupling. Our analysis of this spatial exchange also 
remaiHS-|Valid in systems that use additional quantum 
dots|L3E3 with universal quantum operations in mind, as 
long as each operation is applied to only two dots at a 
time. 

While the movement of quantum dots requires the pre- 
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cise control of the confining potential, which may be dif- 
ficult to realize experimentally, our analysis shows that 
the construction of such a system is viable. With re- 
alistic parameter values from current experiments, our 
analysis produces infidelities smaller than 10~ 4 after ten 
decoupling cycles. This setup could also be a productive 
step towards the experimental realization of more com- 
plicated exchange systems, with many more interesting 



applications. 
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